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A  BOUNDARY  INTEGRAL  METHOD  FOR  EDDY  CURRENT  FLOW  AROUND  CRACKS  IN  THIN  PLATES 


M.A.  Morjeria,  S.  Mukherjee  and  F.C.  Moon 
Department  of  Theoretical  and  Appllad  Mechanics 
Cornell  University,  Ithaca,  New  York 


A  boundary  element  method  which  employs  a  Green's  function  for  a  crack  has  been  developed  to  calculate  the 
induced  eddy  current  flow  around  cracks  in  thin  conducting  placee.  The  theoretical  equations  employ  a  etreaa 
function  for  the  current  density  vector  and  is  equivalent  to  cha  electric  field  vector  potential  method.  A  low 
frequency  or  large  skin  depth  approximation  leads  to  e  Poisson  equation  for  steady  harmonic  inductor  fields. 
Induced  currents  eround  a  creek  in  a  square  place  due  co  a  uniform  inductor  field  for  various  crack  positions 
and  eices  have  been  calculated  in  this  paper. 

The  effect  of  the  relative  position  and  length  of  the  crack,  with  reepecc  to  the  plate  width,  on  the  eddy 
current  density  near  the  cipe  of  the  crack  is  given  special  attention.  These  results  may  be  useful  co  simulate 
eddy  current  flow  detection  phenomena. 

INTRODUCTION 


The  boundary  element  machod  (BEM)  (also  called 
the  boundary  integral  equation  method)  has  emerged  as 
an  important  computational  technique  for  alectrodynamic 
problems.  Wu  at  al  (1]  and  Aacalla  at  al  (2]  have  ad¬ 
dressed  magnetostatic  problems  by  the  BEM  while  Trov~- 
bridge  [31  has  considered  problems  by  cha  magnetic  po¬ 
tential  method.  Very  recently,  Salon  and  Schneider 
[41  have  solved  problems  of  eddy  current  flow  in  long 
prismatic  conductors  by  the  BtM  based  on  Che  electric 
potential  approach. 

In  this  paper,  we  describe  a  powerful  boundary 
element  technique  for  calculating  induced  eddy  currant 
flows  in  conducting  pieces  with  through  cracks  using 
the  electric  potential  approach.  The  BOt  has  the  im¬ 
portant  advantage  that  only  the  boundary  of  a  body 
(.rather  than  the  entire  domain)  needs  to  be  discrecised 
in  a  numerical  solution  procedure. 

There  have  been  soma  attempts  co  model  eddy  cur¬ 
rent  flow  around  annular  cracks  in  rods  and  in  places 
by  replacing  cracks  by  slots  (eee  for  example  Ref. 

[3]).  However,  we  have  shown  thee  the  induced  current 
in  che  vicinity  of  a  crack  leads  to  a  singularity  of 
current  dene icy  at  che  crack  tipe  [6,7].  This  high 
concentration  allows  one  to  use  eddy  current  testing 
devices  such  as  active  and  passive  search  coils  to 
detect  che  presence  of  cracks.  It  also  results  in  a 
temperature  hot  spot  which  can  be  detected  by  infrared 
scanning  [6,3],  The  boundary  element  technique  intro¬ 
duced  by  che  authors  [6,7]  and  described  here  allows 
one  co  model  exactly  che  singular  nature  of  current 


density  at  crack  tips  of  thin  places.  This  technique 
can  handle  any  arbitrary  shape  of  the  plate  and  gener¬ 
al  magnetic  fields. 

In  this  paper  we  discuss  application  of  che  BEM 
to  eddy  current  flow  in  e  cracked  square  plate  due  to 
an  uniform  inductor  field  applied  normal  to  the  plate. 
A  number  of  crack  sizes  Co  plate  size  configurations 
hae  bean  considered.  Also,  effect  of  the  relative 
position  of  a  crack  tip  Co  the  plate  edge  on  the  in¬ 
duced  eddy  current  dietribution  has  been  investigated. 

GOVERNING  EQUATIONS 

A  thin  piece  with  a -crack  in  it  is  shown  in  Fig. 
1.  The  plate  is  made  of  e  conducting  materiel  of 
conductivity  o.  The  plate  boundary  can  be  arbitrary 
and  its  thickness  (uniform) -is  h.  The  thin  line 
crack  is  of  length  2e  end  can  have  arbitrary  ori¬ 
entation  relative  co  the  boundary  of  the  place.  The 
coordinate  system  for  che  problem  is  also  shown  in 
Fig.  1.  The  origin  of  coordinates  lies  at  che  center 
of  che  crack  and  at  che  mideurface  o'  the  plate. 

An  external,  oscillatory  magnetic  field,  g°, 
is  applied  which  induces  e  current  density  J  in  the 
place.  It  is  assumed  chat  the  current  density  is 
uniform  across  che  plate  thickness  end  chat  che  skin 
depth  (vhich  le  Inversely  proportional  co  the  square 
root  of  che  frequency)  is  large  compared  to  che  piece 
thickness. 


(7) 


A  (Cream  function  (or  eleccrlc  potential)  formu¬ 
lation  la  uaed  in  chls  problem.  The  stream  function, 
dx^  ,x,) ,  is  defined  as 

J  -  7*(*k)  .  -Ic  «  7*  (1) 

This  equation  guarantees  the  conservation  of  charge 
equation  7*J  •  0  for  charge  free  regions. 

Using  Ohm's  law  che  governing  differencial  equa¬ 
tion  for  the  stream  function  is  obtained  as  [6,7] 

-  r  £  (B$«J)  <2) 

In  Che  above,  B^  is  Che  self  magnetic  field 
due  to  the  current  J.  It  has  been  shown  in  ref. 

[9],  however,  that  for  a  sinusoidal  applied  field, 
with  che  skin  depth  ouch  greacer  chan  the  thickness 
of  the  plate,  3*  can  be  neglecced  relative  to  the 
applied  field  3°.  This  assumption  simplifies  che 
problem,  and,  with  3^  *  B^  ei“C  (with  i  •  ^-1 
and  -j  the  frequency) ,  the  spatial  pert  of  j,  satis¬ 
fies  a  two-dimensional  nonhomogcncous  Polsacn's 
equation 

7Zp  *  i^oBj  -  f(xltx2)  (3) 

The  boundary  condition  requires  chat  the  current 
must  be  tangential  to  che  plate  boundary.  Thus  4 
is  required  to  be  conscanc  on  che  boundaries  jC^ 
and  jC2 .  On  one  boundary,  the  value  of  9  la  set 
to  zero,  while  on  che  other  boundary  j,  •  C  and  C 
is  obtained  from  che  assumption  chat  che  nec  flux  flow¬ 
ing  through  the  crack  boundary  is  zero.  This  leads  to 
the  condition 

<J  J.tda  •  0  (4) 

3tl  *  ’ 

where  c  is  an  unit  tangent  to  and  s  is  the 

distance  measured  along  e  boundary  in  the  anticlock¬ 
wise  sense.  This  formulation  assumes  chat  no  current 
flows  across  the  creek  or  creek  dp  and  leads  to  a 
singularity  of  the  J  field  at  a  crack  clp.  This  is 
analogous  to  the  stress  singularity  in  fracture  me¬ 
chanics.  It  is  possible  that  soma  leakage  of  currenc 
occurs  across  a  crack  clp  and  thus  relieves  the  singu¬ 
larity  in  actual  conductors.  Possible  leakage  of 
currenc  is  noc  considered  in  chls  oeper.  (It  is  noted 
hare  chat  Infrared  scans  of  eddy  current  flow  around 
cracks  do  indeed  show  a  large  increase  in  temperature 
at  the  creek  tips,  indicating  high  current  density 
at  the  creek  tips  [61.) 

In  susmery,  the  boundary  conditions  on  D ,  used 
in  chls  formulation,  are 

*  *  0  on  the  crack  boundary  3C^  (5) 

•  0  on  the  outside  boundary  SC,  (6) 


SC. 


dll  , 

d7  ds 


0 


These  boundary  conditions,  together  with  the 
field  equation  (3) ,  constitute  a  well  posed  sroblem. 


BOUNDARY  ELEMENT  FORMULATION 


Integral  equations 

An  integral  equation  formulecion  for  Poisson's 
aquation  (3)  can  be  wrlcten  as  (Fig.  1)  [6,7] 

2r*,(p)  •  j  5k(o,0)G(0)dsn  *  /  K(p,q)f(q)dA  (8) 

SC,  (A 

This  is  a  single  layer  pocendel  formulation 
where  G,  a  source  scrength  function  on  the  oucslde 
boundary,  oust  be  determined  from  che  boundary  condi¬ 
tion  on  it  (equation  9).  The  points  o  (or  P)  and 
q  (or  Q)  are  source  end  field  points,  resoecclvely, 
with  capital  letters  denoting  polncs  on  che  boundary 
of  the  body  end  lower  case  leccers  denoting  polncs 
inside  che  body.  The  area  of  the  body  B  is  denoted 
by  A. 

It  has  been  shown  [6]  chat  y  from  eouaclon  (8) 
with  the  following  kernel  satisfies  che  boundary  con¬ 
ditions  (5)  and  (7)  implicitly. 


K(p,q)  -  Re[p(z,7,zo] 

(9) 

where 

*(z.7.zo)  •  tnd-rj/O  -  tnd-rjT) 

z  i 

ri  *  ~  f-2-  •  >ilil 

(10) 

*  /r 7 

Z.Z^±  ,  |«|  <1 

Re  denotes  the  real  pert  of  che  complex  argument, 

z  and  z  ere  the  source  end  field  oolne  coordinates, 
0 

respectively,  in  complex  notation  end  a  suoeroosed 
bar  denotes,  as  usual,  the  complex  conjugate  of  e  Com¬ 
dex  quantity. 

The  remaining  boundary  condition  (6)  on  the  out¬ 
side  surface  is  satisfied  by  using  e  differentiated 
version  of  (8)  and  caking  the  limit  as  p  inside  B 
approaches  a  ooint  P  on  SCj.  Defining 

H.  -  Im(l±  -  |£)  ,  a,  -  -ReOfJ  *  %  (11) 

1  Sz  *  :Z  )z 

che  boundary  condition  (6)  becomes 

0  -  6  H1(P,Q)n1(P)G(Q)dsQ  *  /AHi(P,q)nl(P)£ (q)dAq(12> 

where  nj^  ere  the  components  of  the  unit  outward  norm¬ 
al  to  SC2  at  some  locally  smooth  point  on  ic. 

The  currenc,  J,  at  e  point  inside  the  bodv  Is 
obtained  from  equations  (1)  and  (8). 


Discretization  of  aquations  end  solution  strategy 

The  outer  boundary  of  the  body,  SC,,  is  divided 
into  Sj  straight  boundary  elements  using  (N^  « 
N2)  boundary  nodes  end  che  interior  of  the  body,  A, 
is  divided  into  n^  triangular  internal  elements. 

A  discretized  version  of  equation  (12)  is 

+  r«1/flA1Hi(V')“i(pM>f('»dA,  (u) 

where  is  che  point  P  where  it  coincides  with 
a  node  M  et  a  center  of  e  boundery  segment  on  3C2 
and  qs^  end  AA^  are  boundary  and  Internal  elements 
respectively. 

A  simple  numerical  scheme  is  used  in  which  the 
source  strengths  G  ere  assumed  to  be  piecewise  uni¬ 
form  on  each  boundary  segment  with  their  values  to  be 
determined  at  the  nodes  which  lie  et  the  cencers  of 
each  segment.  Substitution  of  the  piecewise  uniform 
source  strengths  into  equation  (13)  and  carrying  out 
of  the  necessary  integrations,  analytically  and  num¬ 
erically,  leads  to  an  algebraic  system  of  the  type 

(0)  -  [A] (G)  *  (d)  (14) 

The  coefficients  of  che  matrix  [A]  contain 
boundery  integrals  of  the  kernel.  The  vector  (d) 
contains  contributions  from  che  arae  Integrals  end 
the  vector  (G)  che  unknown  source  strengths  et  the 
boundary  nodes.  The  dimension  of  (G)  depends  only 
on  Che  number  of  boundary  elements  on  3Cj  and  the 
internal  discretization  is  necessary  only  for  the 
evaluation  of  integrals  with  known  integrands. 

The  solution  strategy  is  as  follows.  The  matrix 
[A]  end  vector  (d)  in  equation  (14)  are  first  eval¬ 
uated  by  using  the  appropriate  expressions  for  the 
kernels  and  che  prescribed  function  f  in  equation 
(3).  Equation  (14)  is  solved  for  the  vector  (G). 

This  value  of  (G)  is  now  used  in  e  discretized  ver¬ 
sion  of  equation  (8)  to  obtain  the  values  of  the 
stream  function  *  et  any  point  p.  Finally,  the 
current  vector  et  any  point  Is  obtained  from  aquations 
analogous  to  (8) . 


NUMERICAL  RESULTS 

In  the  numerical  computations,  3°  in  Eq.  (13) 
is  assumed  to  be  e  constant.  Eq.  (3)  can  be  non- 
dlmenalonalized  to  the  form 

72i(x1,x2)  -I  ,  x,  ■  x1/s  (15) 

vhers 


14irB°R 


R 


and  the  skin  depth 


J 


12irB°R 


For  che  results  in  this  paoer  a  -  2.  A  typical 
mesh  for  che  results  for  example  shown  in  Fig.  2d 
has  48  boundary  segments  uniformly  distributed  along 
the  upper  half  (due  to  symmetry)  of  the  boundary 
of  the  plate.  In  order  to  evaluate  the  known  area 
integral  in  Equation  13,  the  internal  area  quadra¬ 
ture  was  used.  It  took  about  300  c.p.u.  secs  on 
IBM  370/168  to  obtain  che  results  in  Fig.  2d. 

The  equation  (15)  is  identical  to  one  relating 
to  the  torsion  of  shafts.  The  BEM  was  verified  by 
comparing  the  numerical  results  for  the  solution  of 
(15)  in  e  square  plate  without  e  crack  to  known 
analytical  results  for  che  torsion  of  a  shaft.  The 
BEM  method  has  also  been  checked  against  e  finite 
element  technique  developed  for  eddy  current  prob¬ 
lems  [10). 


Eddy  current  stream  lines  (p  lines)  are  shown 
in  Figs.  2  and  3  for  e  square  plate  with  e  crack  in 
it.  Fig.  2  (s)  -  (c)  shows  how  the  stream  lines 
are  affected  by  varying  che  size  of  the  plate  while 
keeping  the  crack  size  asms.  Due  to  symmetry  only 
the  upper  half  of  the  place  is  shown  in  Fig.  2.  Fig. 
2  (d)  shows  che  effect  of  moving  the  crack  towards 
one  of  the  plate  edges.  Fig.  3  shows  e  close  up  of 
the  stream  lines  near  right  crack  tip  for  Fig.  2  (c). 
The  crowding  of  stream  lines  near  creek  dps  leads 
to  large  gradient  of  *  end  therefore  large  induced 
currents  in  this  region.  The  local  temperature  is 
proportional  to  the  square  of  che  current  density 
(j.J).  Figure  4  shows  calculated  temperature  scens 
elong  e  line  slightly  above  the  crack  (x2  *  .0125) 
for  the  results  shown  in  Fig.  2.  From  Figs.  4  (e)  - 
(c)  one  can  conclude  chat  as  che  crack  size  Increases 
relative  to  the  plate  size  the  hot  spots  et  crack 
tips  are  more  significant  compered  to  those  et  che 
edges.  The  effect  of  moving  the  creek  near  che  place 
edge  gives  rise  to  significant  hot  spots  as  shown 
in  Fig.  4  (d)  end  (c).  This  becomes  more  apparent 
when  we  look  et  the  'Eddy  Current  Intensity  Factor' 
defined  below.  It  has  been  shown  [6,7]  chat  the 
eddy  current  density  squared  is  inversely  proportion¬ 
al  to  the  distance  r  from  a  creek  tip.  He  can  de¬ 
fine  an  eddy  currant  intensity  factor,  M  as 


Table  1  shows  Che  calculated  values  of  M_T^  for  the 
two  crack  tips  for  the  results  shown  In  Fig.  2. 

It  Is  seen  that  the  value  of  M _  remains  practical¬ 

ly  constant  for  varying  place  sizes.  However  It 
changes  significantly  as  a  crack  dp  Is  brought  near 
an  edge  of  the  plate. 
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